This study intends to present a biomechanical model of the lumbar spine applying the finite ELEMENT method in order to evaluate the behavior of the spine with disc degeneration. The high rates of patients suffering from this phenomenon encouraged us to study the effects of disc degeneration on spinal response. In the proposed method, the entire lumbar spine, including the vertebrae L1 to S1, were simulated. Degeneration of the disc was also modeled in three different ways, decreasing disc height, changing the mechanical properties of the nucleus, and changing the properties of ligaments and collagen fibers. This degeneration was considered simultaneously for both L4-5 and L5-S1 discs, which is referred to as double-LEVEL degeneration in this study. After modeling and applying synthetic loading (bending moments with a follower load), the analysis was performed via ABAQUS software. The results, including intradiscal pressures and the intervertebral rotation, were also compared with experimental data for further verification. The findings of this study illustrate that double-LEVEL disc degeneration reduces intradiscal pressures in L4-5 and L5-S1 discs. However, the intradiscal pressure of a degenerated disc does not change the intradiscal pressure of other adjacent discs. Moreover, in extension and axial rotation loading, increasing disc degeneration would lead to an increase in intervertebral motion.